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Abstract. In second-generation, ground-based interferometric gravitational- wave 
detectors such as Advanced LIGO, the dominant noise at frequencies / ~ 40 Hz to 
~ 200 Hz is expected to be due to thermal fluctuations in the mirrors' substrates and 
coatings which induce random fluctuations in the shape of the mirror face. The laser- 
light beam averages over these fluctuations; the larger the beam and the flatter its 
^ ■ light-power distribution, the better the averaging and the lower the resulting thermal 

noise. In semi-infinite mirrors, scaling laws for the influence of beam shape on the four 
dominant types of thermal noise (coating Brownian, coating thcrmoclastic, substrate 
. Brownian, and substrate thcrmoclastic) have been suggested by various researchers 

and derived with varying degrees of rigour. Because these scaling laws are important 
tools for current research on optimizing the beam shape, it is important to firm up 
our understanding of them. This paper (1) gives a summary of the prior work and of 
gaps in the prior analyses, (2) gives a unified and rigorous derivation of all four scaling 
laws, and (3) explores, relying on work by J. Agrcsti, deviations from the scaling laws 
due to finite mirror size. 
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PACS numbers: 04.80.Nn 
1. Introduction and Summary 

Second generation interferometric gravitational wave detectors such as Advanced LIGO 
will be approximately ten times more sensitive than the current LIGO interferometers, 
leading to an improvement in event rate such that the first few hours of data at design 
sensitivity will contain more signals than the entire year-long science run that is presently 
under way [1]. In advanced LIGO's most sensitive frequency band (/ ~ 40 to 200 Hz), 
the sensitivity is limited by internal thermal noise, i.e., by noise in the substrates and 
reflective coatings of the four test masses (see, e.g., figure 1 of [2]). Lowering the internal 
thermal noise would increase advanced LIGO's event rate throughout that band. 

Internal thermal noise can be divided into two different types: Brownian thermal 
noise (due to imperfections in the substrate or coating material, which couple normal 
modes of vibration to each other) and thermoelastic noise (due to random flow of heat 
in the substrate or coating, which causes random thermal expansion). When the laser 
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beam shape is Gaussian, the Brownian and thermoelastic noises in the substrate (e.g. 
[3]) and in the coating (e.g. [I] [5]) are well understood. 

One way of lowering the internal thermal noise is to i) flatten the shape of the laser 
beam that measures the test mass position so it better averages over the mirror faces' 
fluctuating shapes, and ii) enlarge it to the largest size permitted by diffraction losses. A 
specific enlarged, flattened shape, the mesa beam, has been proposed by O'Shaughnessy 
and Thorne and explored (theoretically) in detail by them, d'Ambrosio, Strigin and 
Vyatchanin [6j [7J E], and by Agresti and DeSalvo 0, HO]. The mesa shape was found to 
reduce the thermal noise powers by factors of order two, with corresponding significant 
increases in the distances to which the planned interferometers can search. Motivated 
by this, mesa beams are currently being explored experimentally |llj|12j . 

The mesa shape is unlikely to be optimal. Bondarescu and Chen (Caltech/AEI) 
are currently seeking the optimal beam shape for each of the four types of noise; they 
are also seeking a balance between the competing demands of the four optimal shapes. 
Further research will require balancing practical aspects of mirror design against the 
(possibly impractical) ideal shapes. 

In all this current research, a crucial tool is a set of scaling laws for the dependence 
of the four types of thermal noise on the beam shape, in the limit of a mirror that is 
large compared to the beam diameter ( "semi- infinite mirror"). These scaling laws have 
been proposed by various researchers over the past several years, and they have been 
derived with varying degrees of rigour, and in some cases with unnecessarily restrictive 
assumptions. This prior work will be discussed and critiqued in section 12.21 

Because these scaling laws are so important for current research, this paper 
scrutinizes them and their accuracies in some detail. In section [2] the scaling laws and 
assumptions underlying them are presented and prior research on them is described. 
Then in section [3] a unified and rigorous derivation of all four scaling laws is presented. 
In section H] the breakdown of the scaling laws due to finite mirror size is explored. And 
finally, in section [5] a few conclusions are given. 

2. The scaling laws and prior research on them 

2.1. Model and Summary 

To explore the effect of the beam shape on the internal thermal noise, I consider a 
cylindrical test mass substrate of radius R and thickness H and suppose that these size 
scales are comparable: R ~ H. I choose a cylindrical coordinate system (r, tp, z) such 
that r = is the mirror axis, z = is the reflectively coated surface of the mirror 
substrate, and points with < z < H are inside the mirror substrate. 

An axisymmetric laser beam with intensity profile p(r) is normally incident on the 



Dependence of test-mass thermal noises on beam shape in GW interferometers 



mirroiyj. The intensity profile is normalized, so 
f R 

2ir / drrp{r) = 1. (1) 
Jo 

The beam measures q(t), a weighted average of the mirror's longitudinal position 
Z(r,<f,t) (equation (3) of [T3] ) 

r2n i-B, 

q(t) = / dip / drrp(r)Z(r, ip,t). (2) 
Jo Jo 

In LIGO, so as to keep diffraction losses < lppm, the beam radius over which, say, 
95% of the signal q(t) is collected, is kept significantly smaller than the mirror radius R 
and thickness H. This motivates the idealization of the mirror as a semi-infinite slab 
bounded by a plane, R — > oo, H — > oo. (The accuracy of this infinite-test-mass (ITM) 
approximation will be discussed in section I4~2l ) 

Internal thermal noise will cause small fluctuations in the longitudinal position of 
the mirror Z(r,(p,t). The spectral density S q associated with the measurement of the 
mirror position q is given by the fluctuation dissipation theorem (equation (1) of [13]): 

_ 2ksT Wdiss 

b « ~ n 2 p F 2 ■ W 

Here ks is Boltzmann's constant, T is the temperature of the material, and Wdiss is the 
power that would be dissipated if a longitudinal force F with frequency / and pressure 
distribution p(r) were applied to the mirror surface (Levin's [13] thought experiment). 
Because the frequencies of interest (i.e. / ~ 100 Hz) are far below the lowest resonant 
frequencies of the mirror / res ~ (a few km/s)/ (about 10 cm) ~ 10 Hz, the hypothetical 
applied force F can be idealized as static when computing the resulting strain of the 
mirror. 

Thus the noise S q can be computed using the following algorithm: 

(i) Statically deform the (semi-infinite) mirror with a force F with pressure distribution 
p(r) the same as the light's intensity profile; 

(ii) compute the Brownian and thermoelastic dissipated power Wdiss due to the 
deformation caused by F; 

(iii) substitute H^iss into equation ([3]) to get the spectral density S q of the thermal noise 
of a measurement of the average position q. 

Note that from S q , one can easily compute the thermally-induced gravitational- 
wave-strain noise power Sh{f) in a measurement by the interferometer. If mirrors 1 and 
2 are in one arm (of length L = 4km), and mirrors 3 and 4 are in the other arm (also 
of length L), the interferometer measures h = [(qi — ^2) — (93 — q±)]/L, where 91,2,3,4 are 
the measured positions of the four mirrors. Because the noises in the four test masses 
are uncorrelated, the spectral density is just = (4/L 2 )S q . In the remainder of 

I The shape of the mirror faces must also be changed slightly (by height changes < one wavelength of 
the laser light) so that p(r) is an cigenmodc of the arm cavity. In this paper, I assume that the mirror 
faces take whatever shape is necessary to support a beam with intensity profile p(r). 
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this article, when referring to the noise of a single test mass, the subscript "q" will be 
suppressed (i.e. S = S q ), while the gravitational- wave-strain noise power will always be 
referred to as Sh- 
in section |3~T| I compute the strain distribution that results from applying the force 
F to a homogeneous, isotropic, semi-infinite mirror with a very thin reflective coating 
of a possibly different material. The calculation is a straightforward generalization of 
section 2 of [1] . In this calculation, I model the coating as a thin layer (of order microns, 
as compared to the cm size scales of the substrate) which adheres to the mirror surface. 
In section 13.21 1 use the strain distributions to compute each of the four types of thermal 
noise S(f). I find that if pi(r) and p 2 {r) are two different beam shapes, then 



where n = 1 for coating Brownian and coating thermoelastic noise, n = for substrate 
Brownian noise, and n = 2 for substrate thermoelastic noise. Here p(k) is (up to factors 
of 2n) the two-dimensional Fourier transform of p(r) over the surface of the mirror: 



Here Jq(x) is the th Bessel function of the first kind (the axisymmetry allows the 2D 
Fourier transform to reduce to a ID Hankel transform). 

If one knows Si )Tl , computing Sa.n amounts to computing simple integrals of p(k). 
If one holds everything else fixed but changes the beam shape, the scaling law (J4]) 
makes it straightforward to determine the improvement in the thermal noises and the 
corresponding improvement in the interferometer sensitivity. 

In the remainder of this paper, I derive these scaling laws, comment on their 
implications for advanced LIGO, and estimate their accuracy for finite test masses. 
In section 12.21 I discuss prior work related to the scaling laws. In section 13.11 I 
compute the strain SV, due to a hypothetical applied force F with pressure distribution 
p(r). Then, in section I3T21 I compute the dissipated power Wdi SS for the Brownian and 
thermoelastic dissipation in the coating and the substrate and insert Wdiss i n t° equation 
([3]) to determine how the noise depends on the beam shape. In section 14.11 I discuss 
implications of this result for advanced LIGO, and in section |4~21 1 discuss the accuracies 
of the infinite-test-mass (ITM) scaling laws by comparing with others' finite-test-mass 
(FTM) predictions for the cases of Gaussian and mesa beam shapes. I make some 
concluding remarks in section [5l 

2.2. Discussion of prior research 



Si,« _ JTdM-^)] 2 



(4) 



S 2 ,n j^dkk n [p 2 (k)} 2 




(5) 



2.2.1. Thermoelastic substrate noise In Levin's thought experiment, the dissipation 
associated with thermoelastic noise arises from heat flow down temperature gradients, 
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which are induced by compression of the coating or substrate by the force F. The 
increase in entropy corresponds to a dissipated power. 

In 2003, in connection with his invention of the mesa beam and exploration of 
its properties, O'Shaughnessy used Levin's thought experiment to derive the following 
scaling law for the thermoelastic substrate noise: 



where does not depend on the beam shape. This scaling law ultimately motivated 
the other three. O'Shaughnessy included his derivation of this law as Appendix H of his 
(as yet unpublished) 2004 paper with Strigin and Vyatchanin [8] on mesa beams. He 
used a slightly different (but no less rigorous) method from the unified derivation I give 
in section 13.2.41 [equation (1381) ] . O'Shaughnessy wrote the scaling law in terms of 2D 
Fourier transforms; the reduction to ID Hankel transforms makes numerical evaluation 
of the scaling law (El) very efficient (section [4.1]) . 

2.2.2. Thermoelastic coating noise Braginsky and Vyatchanin (Appendix B.2 of [5]) 
and Fejer and collaborators (section IV D of [H]) have independently calculated the 
thermoelastic coating noise for Gaussian beam shapes (though the analysis in [5] is only 
valid when the coating and substrate elastic properties are identical [E]). Scrutinizing 
the derivation in [H], Thorne speculated in 2004 (unpublished) that the thermoelastic 
coating noise obeys a scaling law of the form 



In 2006 I verified Thome's conjecture via almost trivial generalizations of the Braginsky- 
Vyatchanin and Fejer et. al. analyses; my derivation is given in section 13.2.21 [equation 
( [28]) ]. In 2006 O'Shaughnessy, learning of my work but not knowing my result, extended 
a clever dimensional analysis argument (section 12.2.51) that he originally invented 
for Brownian coating noise (below) to the other three types of noises pj)]; but for 
thermoelastic coating noise he got an answer that disagrees with Thome's conjecture 
and my result (ED). When I pointed out the discrepancy, O'Shaughnessy found an error 
in his dimensional analysis and revised his manuscript to give the correct scaling law 



2.2.3. Brownian coating noise In Levin's thought experiment, the dissipation 
associated with Brownian thermal noise can be modelled as arising from a loss angle, 
which is an imaginary (i.e. damping) correction to the material's Young's modulus 
caused by coating or substrate imperfections. 

In 2004, Thorne communicated to O'Shaughnessy and Vyatchanin his conjecture 
([7]) for the scaling law for thermoelastic coating noise, and challenged them to find an 
analogous scaling law for Brownian coating noise. Independently, they each devised 




(6) 




(7) 



© [16j. 
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simple arguments that led to the law 




coat 



/ dkkp\k) 



(8) 



O'Shaughnessy gave both an argument based on dimensional analysis (section [2.2.5p and 
a derivation for the special case that the substrate and coating have the same elastic 
properties. Vyatchanin's analysis pT] was based on a derivation for Gaussian beams, 
followed by an argument that, if a beam with another shape p(r) can be constructed by 
superposing Gaussian beams, then this scaling law must hold also for that other shape. 

The scaling law (jSj) is local, i.e., the noise at a point on the mirror depends only 
on the beam intensity evaluated at that point. Thome's intuition, however, led him 
to believe (incorrectly) that the scaling law should be non/ocaf§|. Consequently, Thorne 
was so highly sceptical of O'Shaughnessy's and Vyatchanin's arguments and the claimed 
scaling law that he — unfortunately — dissuaded both O'Shaughnessy and Vyatchanin 
from publishing their arguments and result. 

The following year (2005), Thorne, still sceptical of the O'Shaughnessy- Vyatchanin 
result (JSj), suggested to me that I carry out a detailed derivation of the Brownian- 
coating-noise scaling law from first principles. My analysis, based on Levin's method 
and reported in this paper, gave the result (jHJ), in agreement with O'Shaughnessy 
and Vyatchanin, and motivated O'Shaughnessy to publish [16] his dimensional-analysis 
argument. 

O'Shaughnessy's derivation is restricted (unrealistically) to identical elastic 
properties for substrate and coating. My derivation [equation ( |23l) below] allows the 
substrate and the coating to have different elastic properties. Vyatchanin's derivation 
is valid only for those beam shapes that can be achieved by superposing Gaussians - 
though it might well be that any shape can be achieved in this way. My derivation is 
definitely valid for any axially symmetric beam shape p(r). 

2.2.4- Brownian substrate noise In 2005 Vinet proposed [18] the following scaling law 
for the substrate Brownian noise: 



He deduced this law as a trivial consequence of his equations (1) - (3). He did not 
present a derivation of those equations, but he recognized that they can be obtained by 
generalizing the derivation in [19], which assumes that the beam shape is Gaussian. In 
section [3.2.31 I explicitly derive equation (Q. In parallel with my work, O'Shaughnessy 
applied his dimensional analysis technique to verify Vinet's scaling law ([91). 

2.2.5. Dimensional analysis O'Shaughnessy's dimensional analysis argument, referred 
to above, consists of three steps: 




(9) 



§ It turns out (section !3.2.ip that nonlocal terms do appear at intermediate steps in the derivation but 
do not contribute to the scaling law itself. 
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(i) The scaling laws must take the form of a translation-invariant inner product of p(r) 
with itself, since the mirror is taken to be semi-infinite. In the Fourier domain, 
for axisymmetric beam shapes pi(r) and P2{r), the scaling law must then take the 
form: 

51 = J °° dkkG(k)pl(k) 

5 2 J °° dkkG(k)p 2 2 (k)' [ ' 

(ii) The only length scale (other than the width of the beam) is the small coating 
thickness d, so G(k) = k n d for coating thermal noise and G(k) = k n for substrate 
thermal noise. 

(iii) The power n is chosen by demanding that, when the beam shape is a Gaussian, the 
noise scale as the correct power of the beam width. 

This argument turns out to produce the correct scaling laws, but without sufficient 
care, it can also lead one amiss. For instance, when considering thermoelastic coating 
noise, step (ii) must be amended, since there is a second length scale: the characteristic 
length of diffusive heat flow [5][T3]. In his original manuscript [15], O'Shaughnessy 
neglected this second length scale, and incorrectly deduced that n = 3 for coating 
thermoelastic noise. After I contacted O'Shaughnessy regarding this error, he corrected 
his analysis [16] and obtained the same result, n = 1, as I had derived (section 13.2.21) 
below. 



3. Derivation of the infinite-test-mass (ITM) scaling laws 

3.1. Strain of a semi-infinite body with thin facial coatings due to a static, 
axisymmetric force 

The thermal noise is determined by the symmetric part of the strain Sij that the test 
mass would experience if a normal force with pressure p(r) were applied to the mirror 
surface. In this section, I evaluate S^ in the mirror substrate and coating. In section I3T21 
I use these results to compute Wdiss [which, by equation ([3]), determines the thermal 
noise] . 

If the displacement vector of an element of the test mass is Ui, then the strain S^ is 
= VjUi. Following the methods developed in [19] (but correcting some typographical 
errors), equation (19) of [3] gives the cylindrical components of the displacement of the 
test mass substrate: 

Ur = — r dkJAkr)e- kz ( 1 - + kz] p(k), 

2^ Jo \ X + fi J 

(11a) 

u 9 = 0, (116) 



u z = ±- dkJ (kr)e- kz (l + + kz) p(k). 



:uc) 
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Here A and \x are the Lame coefficients of the substrate. The vector Ui satisfies the 
equilibrium equation VjTjj = 0. (Throughout this paper, I use the Einstein summation 
convention.) 

The non- vanishing components of the symmetric part of the strain are [with commas 
denoting partial derivatives] 

9 = Su, (12a) 
S rr = u r ^ r = 9 S zz Sptp, (126) 

= ~i ( 12c ) 

S(rz) = S( zr ) = 7^{u ryZ + U 2 , r ), (12d) 

S zz = u z , z . (12e) 

Evaluating the derivatives of equations filial) - (JTTcj) and inserting the result into 
equations (jl2aj) - (11 2 el) gives 



9 = -L f°° dkkJ (kr) (yT~) e- kz P(k), 
2/Wo V A + AV 



Srr 9 S zz Stpip) 



kz p(k), 



S„ = 1 fd^Wl-^ 

W 2/i7 r V ^ + M 

S {zr) = - — / dkkJ^kr^kz^pik), 

2yU Jo 

^ = ^ r dkkJ {kr) (--^— - kz) e~ kz p(k). 



(13a) 
(136) 

(13c) 

(13a 1 ) 

(13 e) 



Setting z = in equations (jl3ap - (jl3e|) and combining with equation (JSJ) yields the 
nonvanishing stresses on the substrate surface: 

^rr | 2=0 = J f TT") ^ ~ S ^ l*=0 ' ( 146 ) 

* w U = J ( T7: ) r d*^2p(*), (14c) 



A - 


- pL 




1 


A- 


- pL 




1 


A- 


" Z 1 



^L=o =o( — IpW- 

Here I have used the identity 

dkkJ (kr)J {kr') = ~ T \ (15) 

r 
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Note that on the surface of the substrate 9 and S zz are local [i.e. their values at any point 
depend only on the value of p(r) at that point], while S w is nonlocal. The component 
S rr can be written as the sum of a local part and a nonlocal part; the nonlocal part of 
S rr is just S^. 

The thin coating approximation gives the nonvanishing components of the coating 
strain in terms of the strain on the substrate surface (equation (A4) of 0]): 

A + 2// coat 



6 



coat 



Acoat 2/X coa t 

A + 2/i 



S 



zz)\ z =0 



Acoat + 



SCO 



coat 



s, 



rr\ z =0 



2/^coat 
acoat _ 



,5' 



zz\z=0 ' 



ocoat ocoat 

ipip ' zz 1 



s. 



coat 

tptp 



s, 



s 



coat 

zz 



Wlz=0 ' 



A - A r 



16a) 



;i6t) 



'16c) 



Acoat 2/l CO at 

A + 2/i 



-S 



Acoat + 2/X c 



S Z z 



zz)\z=0 



z=0 ' 



(i6d) 

In jl] , these conditions are said to hold in the limit that the Poisson ratio of the substrate 



and coating are "not too different," but this restriction is unnecessary (see Appendix B ) . 

Finally, after inserting equations (114q|) - (114 a]) into equations (jlfiqp - (\16d\i I 
conclude that # coat and S z ° at are local; while S^f* and S^° at are nonlocal. However, 
this nonlocality turns out not to influence the coating noises. This is because, after 
using equation ()166|) to eliminate S*™ at , it turns out that the remaining nonlocal part 
gcoat on jy a pp ears j n the coating H^diss [according to equations fl22|) and (127)) ] via the 
integral 



A rr ocoat ocoat 



drr 



( ocoat \ 2 I / ocoat \ 2 i / ocoat \ 
{O rr ) -f- ) ■+■ {S zz ) 



drr 



//icoat ocoat \ ^ i / ocoat \ 2 i of C 

[U - b zz ) + [b zz ) + 2{b, 



coat \ 2 
(pip I 



2^ocoat ^coat 



ST*) 



;i7) 



In Appendix A I show that 



drr(S™ at ) 2 - S™ at (# coat - S™ at ) = 0, (18) 

so only the local parts of the strain (# coat and 5'™ at ) influence the thermal noise. This 
fact turns out to imply local coating scaling laws in agreement with O'Shaughnessy's 
[T6] and Vyatchanin's [UJ arguments (section [372]) . 



3.2. Internal thermal noise 
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3.2.1. Brownian coating noise For Brownian thermal noise in an elastic material, 
the dissipated power is [equation (12) of [13] with a static applied force and with 
U = -(1/2)^] 

rd p2tv poo 

W diss = -nf dz dip drr0(/)<%% (19) 
Jo Jo Jo 

Here <fi is the loss angle (i.e., the imaginary, damping part of the Young's modulus of the 
coating material) and Tij is the stress. When the material is the thin reflective coating 
of a mirror, there are effectively two loss angles [3], <f>\\ and <j)±, defined so that in the 
previous equation 

= (f)s^T^ + (<j> ± - 0n) stt:T- (20) 

This result can be obtained by combining equations (4) and (13) - (15) of [4] with 
equation (9) of [13] and recalling that in the coating, the strain [equations (1 1 6 <aj) - 
( 116ril) ] is diagonal. 

For a homogeneous coating, the stress T^° at is 

rpco&t ^coat^coat^ 2^j coa t ( Cj' c p a t (21) 

where A coat and /i coat are the Lame coefficients of the coating, S?°^ is the symmetric 
part of the coating strain, and 9 = Sa is the expansion. Combining equations (1201 . ffl9l) 
and ff2T]) gives the following expression for W^ss- 

poo poo 

H/ diss = 27r 2 /^||(/) / drrA + 27T 2 fd[Mf)-H(f)] / drrB , 
Jo Jo 

a ( \coat/i2 I r> coat ocoat ocoat\ 

n — \ A ^coat ' Z t J > ^{ij) ^(ij) ) 5 

jzj ^coat ^^coat^coat _j_ 2^/ coa ^ t S* coa ^) (22) 

Combining equations ( 1221 . (1 16(d) - ( 116dl) . ( 11 4 aft - ( 114dft . and (fTBi) and then inserting 
the result into equation ([3]) gives the spectral density S of the Brownian coating noise. 
However, for the present purpose, only terms involving the beam shape are relevant. 
Absorbing all other terms into a single constant Cf^ t yields 

POO 

C = C drV(r). (23) 
Jo 

This is a local scaling law; i.e., the noise at each point on the mirror's surface is 
proportional to the square of the beam intensity there. This law is the same as 
O'Shaughnessy's [16] and Vyatchanin's [17] scaling law for the Brownian coating thermal 
noise. 

Parseval's equation [which follows from equation ([5])] makes it easy to rewrite this 
scaling law in the Fourier domain, which will facilitate comparison with the substrate 
noise. The result is 

POO 

Sf l = C™ / dkkp 2 (k). (24) 
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3.2.2. Thermoelastic Coating Noise The calculation of the thermoelastic coating noise 
is similar to the calculation of Brownian coating noise. But now, in response to the 
static, normal applied pressure p(r), the dissipated power V^diss is caused by heat flow, 
cx VST, down a temperature gradient VST caused by the material's deformation: 



[equation (5) of [3] in the case of a static applied force and after evaluating the time 
average and trivial <p integral] . Here T is the temperature of the coating in the absence 
of the deformation and k is the material's coefficient of thermal conductivity. 

Braginsky and Vyatchanin [5] and Fejer and collaborators [Hj have independently 
solved for the thermoelastic coating noise. The results obtained in [5] are correct only 
when the coating and substrate have the same elastic properties (section I in |14j); 
however, this restriction is not relevant here, since [HJ and [5] agree on the coating 
thermoelastic noise's dependence on the beam shape p(r). 

If the temperature change were adiabatic, ST would simply be proportional to 6* coat 
(see, e.g., equation (12) of [3]). (Physically, this simply means that the temperature 
of an element in the coating changes linearly with volume.) However, as noted in [5], 
the diffusive heat characteristic length £p of the substrate and coating (on the order of 
mm) is far larger than the coating thickness d (which is on the order of a few microns). 
Because diffusive heat flow in the longitudinal direction is not negligible, heat flow in 
the direction normal to the coating cannot be treated adiabatically [5]. By contrast, 
the substrate thermoelastic noise can be treated adiabatically (section 13. 2. 41) . as can the 
heat flow in the plane of the coating ("tangential" heat flow). 

Because the tangential heat flow is adiabatic, dST/dr ~ 9/w, where w ~ cm is the 
length scale over which p{r) varies. On the other hand, dST/dz ~ 9/£d, where Id ~ mm 
is the diffusive heat characteristic length. Because the tangential derivatives are much 
smaller than the longitudinal derivatives, all derivatives except d/dz may be neglected. 
It follows that Wdiss w iH depend only on p(r) and not on its radial derivatives. 

Based on these observations, Braginsky and Vyatchanin [5] and Fejer and 
collaborators [Hj solve the thermoconductivity equation (e.g., equation (1) of [H]) 
for the temperature perturbations ST. Both [5] and [H] assume that the beam shape 
is Gaussian, but it is quite easy to generalize their arguments to non-Gaussian beam 
shapes. Combining equations (B5) - (B7), (66), and (68) of [H] (but now regarding 
their function p(r) as a generic beam shape) shows that the temperature perturbations 
in the coating have the form 



where F(z) is a function of z only. [Equivalently, equation ( 1261) can be obtained by 
combining equations (B.10) and (B.12) of [5] (but now regarding 9 as an expansion 
corresponding to a generic beam shape) with equations (11 4 aft and (116 a[) .] The precise 
form of F(z) is given in [5] and [H] but is not needed in the present discussion. 




(25) 



ST oc p(r) x F(z) 



(26) 
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Next, Braginsky and Vyatchanin compute the squared gradient (V5T) 2 ~ 
(d5T/dz) 2 in equation fl25l) to obtain Wdissi Fejer and collaborators instead compute 
W'diss by considering the interaction of i) the unperturbed stress and strain [i.e., the 
stress and strain due to p{r) when temperature perturbations are neglected], and ii) 
the (complex) perturbations of the stress and strain caused by the small temperature 
perturbations ST. Both methods lead to the following expression for W^ss- (equations 
(B.13) and (B.10) of [5]; equation (69) of [H]) 

POO 

Wdiss = const x / drrp 2 (r). (27) 
Jo 

Plugging this result into equation ([3]) gives the scaling law 

poo poo 

Sll = CH drrp 2 (r)=Cj l dkkp 2 {k). (28) 
Jo Jo 

This is the same scaling law as for Brownian coating thermal noise. The coating 

thermoelastic noise is local and is the same as O'Shaughnessy's p3] and Vyatchanin's 

[T7] law for Brownian coating thermal noise. 



3.2.3. Brownian Substrate Noise For Brownian substrate thermal noise there is only 
one relevant loss angle, 0, so the dissipated power is (equation (49) of [3] with a static 
applied force) 



POO POO 

W diss = 2tt 2 /0(/) / dz / drr (\6 2 + 2pLS {ij) S (ij) ) 
Jo Jo 



(29) 



The integral of the squared strain can be expanded as 



pOC poo poo poo 

/ dz drrS {ij) S {ij) = dz drr (S 2 r + S 2 W + S 2 ZZ ) 
Jo Jo Jo Jo 

poo poo 

= dz drr[(0- S zz f + S 2 ZZ + 2S 2 rz) + 2S 2 W - 2S W (9 - S„)](3Q) 
Jo Jo 



In Appendix A , I show that 

POO 

/ drr [Si - S w ( 
Jo 

Substituting this result into equation (1301) yields 



S z 



OO POO 



(31) 



(32) 



Wdiss = 2tt 2 /0(/) / dz / drr X9 +2^(6- S zz ) z + 2^S 2 ZZ + A^iS( rz) 
Jo Jo 1 

This expression can be evaluated term by term. Inserting equation (113al) into the integral 
of 9 2 gives 



OO POO 



dz \ drr9 2 
Jo 



1 / 2fi 



V \X + ti 



poo poo poo poo 

I dkkp(k) / dk'k'p{k') / dze~ {k+k ' )z / drrJ (kr)J (k'r). (33) 
o Jo Jo Jo 
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Using the identity 



o 



drrJ n {kr)J n {k'r) = ^— '- (34) 



on equation (|33|) and evaluating the integral over z yields 

J e = J-Jt^-) 2 fWW- ( 35 ) 



o 



ifi 2 \X + \i t 

The other terms in equation ( 1321) can be evaluated similarly; they all turn out to have 
the same dependence on p(k) as Ig. Inserting this result for Wdi SS into equation ([3]) gives 
the scaling law 



r>00 

S?£ = Cifl dkp\k). (36) 



o 



This scaling law is the same as the scaling law ( 1241) for the coating thermal noise except 
that the z integration has reduced the power of k in the integrand by one. This scaling 
law agrees with equations (1) - (2) of [IS]. 



3.2.4- Thermoelastic Substrate Noise In contrast to the case of coating thermoelastic 
noise, the substrate thermoelastic noise can be treated using the adiabatic 
approximation. Therefore, the temperature perturbations 5T that drive the substrate 
thermoelastic noise S"™ are proportional to the expansion, i.e. 5T oc 9. This implies 
(e.g., by equation (j25|) . or equation (13) of [3]) 

roo roo 

SJ1 = Cj u l / d^ / drr(W) 2 (37) 
Jo Jo 

with independent of the strain (and thus also the beam shape). Inserting equation 
(I13ap into equation (137|) gives the scaling law; after absorbing all constants into C™; h 
takes the form 

POD 

SJ2 = CH dkkY(k), (38) 
Jo 

which O'Shaughnessy, Strigin, and Vyatchanin obtain in [8]. This scaling law is the 
same as the scaling law (|36|) for the substrate Brownian noise except that the gradient 
raises the power of k by two. 



4. Applying the ITM scaling laws to second-generation gravitational-wave 
interferometers 

To illustrate the scaling laws (123!) . (j28|) . (136|) . and (1381) . suppose that the noise SV,fc 
[with beam shape pk{r) and thermal noise type r] is known. Here and throughout the 
remainder of this article, r is a label that takes one of the following values: coating 
Brownian (Coat BR), coating thermoelastic (Coat TE), substrate Brownian (Sub BR), 
or substrate thermoelastic (Sub TE). 
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Figure 1. A plot of p ga uss( r /fr) an( i Pmcsa( r /b) for beams with 1 ppm diffraction 
losses (in the clipping approximation) on a mirror with radius R = 17cm. Here 
b = y/ L\/2tt = 2.6cm is the width of the smallest Gaussian beam that can resonate 
in a LIGO arm cavity with length L = 4km and light wavelength A = 1064nm. 



Now, if the beam shape were changed to p u (r) [while holding everything}]] else fixed], 
then the unknown noise S T>U [with beam shape p u {f)] would be [equation (j4j)]: 

S T ,u = C*ITM [T]Pu,Pk]S T ,k, (39) 



with 



and 



ClTM [ T 'iPu,Pk\ 



n(r) 



J™dkk n ^p 2 u (k) 
J °° dkk n (-)p 2 k (k) 

t = CoatBRorCoatTE 
r = SubBR 
r = SubTE 



(40) 



(41) 



When the beam shape is changed from p^ to p u , the amplitude sensitivity changes by a 
factor of Ci TM [r;p u ,Pk]- 



4-1. Implications for advanced LIGO 

In advanced LIGO, the thermal noise may be significantly reduced by changing the 
shape of the laser beam. One proposal is to replace the Gaussian beam shape with a 
mesa beam (also called a flat-top beam) [6]. O'Shaughnessy, Strigin, and Vyatchanin 
[8] have calculated the resulting reduction in substrate thermoelastic noise, Vinet has 
done the same for substrate Brownian thermal noise [18] and Agresti [2] and Agresti 
and DeSalvo [9], [TO] have done the same for both substrate and coating thermal noises — 
all for the realistic case of finite mirrors. The reduction in thermal noise can also be 

|| Since here I am neglecting edge effects, "everything" means the temperature, the materials' elastic 
and thermal properties, the coating thickness, and the frequency. In section l4~2l when edge effects are 
considered, it will be the diffraction loss, not the mirror size, that is held fixed. 
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Figure 2. A log-log plot of the Gaussian beam- width parameter w and mesa beam- 
width parameter D as a function of mirror radius R (top of figure), for mirrors with 
1 ppm diffraction loss in the clipping approximation. The ratio D/w is shown on the 
bottom of the figure. The parameter b is defined in figured] 



understood as a consequence of the simple ITM scaling laws derived above. Although 
I only compare Gaussian and flat-top beams here, the scaling law given in equation 
(JH) makes it simple — if one neglects finite-test-mass (FTM) effects — to compute the 
relative change in sensitivity for any two beam shapes. 
The normalized Gaussian beam shape is 



Pgauss(w;r) 



-r 2 /w 2 



(42) 



where w is the width of the Gaussian beam. It is straightforward to compute p gauS s(^; k), 
since the integral can be done analytically; the result is 



V 



gauss 



[w-k) 



drrJ (kr) 



-r 2 /w 2 



_ c -k 2 w 2 /4 

2tt 



In position space, the mesa beam can be written as (equation (2.5) of 



i{D;r) = N 



i) 



2n I dr r exp 



r 2 + r /2 )(l-i) 
2b 2 



X In 



rr'(l — i) 
b 2 



(43) 



(44) 



Here D is a measure of the width of the beam, b = \J \L/2n, with L = 4km the arm 
length and A = 1064nm the wavelength of the laser beam's primary frequency, and N is 
a normalization constant adjusted so equation (JTJ is satisfied. Note that p meS a( r ) must 
be evaluated numerically; to compute p(k) efficiently, I use the Fast Hankel Transform 
algorithm [20] . 

Examples of the Gaussian and mesa shapes are plotted in figure [TJ In figure the 
width parameters w and D of a sequencd^ of Gaussian and mesa beams are plotted as 
a function of mirror radius R for beams with 1 ppm of diffraction loss in the clipping 

This particular sequence was chosen to facilitate comparison with the results of [10] . which includes 
edge effects. 
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Figure 3. The scaling of thermal noises with beam width w for Gaussian beams 
in the infinite-test-mass (ITM) approximation. More specifically: a log-log plot of 

ClTM [ T ] Pgauss 0;r) 

j Pgauss 

(w a ;r)] as a function of w/b. Here w a /b = 1.24, which 
corresponds to R = 12cm and 1 ppm diffraction losses. Each curve is a power law 
obeying C oc 1/w 1 '. 

approximationEl The ratio D/w is also shown on the bottom horizontal axis. It is 
sometimes useful to regard D and w (for 1 ppm losses) as functions of D/w rather than 
of R — with D/w actually being a surrogate for R. 

The following three cases use equations ( 1391) - ( 1411) to illustrate how the thermal 
noise in advanced LIGO changes with different choices of Gaussian and mesa beam 
shapes. 

4-1.1. Noise of a resized Gaussian beam Suppose puir) = p ga uss(w ; r). Then the 
thermal noises for a Gaussian beam of some different size w are determined by evaluating 
Citm[t; P g auss(w; r), p gauss (ii; ; r)] [equation (14*01) ] and inserting the result into equation 
( 139]) . In this well-known case (see, e.g., the discussion and references in [TO]). Citm can 
be evaluated analytically, yielding the following relation: 

ClTM^; Pgauss(w; r), Pgauss(w ; T*)] OC ■ (45) 

In figure El Citm[t; P g auss(^; r),p ganss (w ] r)] is plotted as a function of the beam width 
w. 

4-1.2. Noise of a resized mesa beam Suppose Pfe(r) = p m esa(-D ; r). Then the 
thermal noises for a mesa beam of some different size D are determined by evaluating 
CiTM[r;Pmesa(-D; r) , p mcsa (_D ; r)] [equation ( 14*01 ] and inserting the results into equation 
("39]) . As shown in figure HI in this case Citm does not scale as an exact power of D 
(although the actual relations are very well approximated by power laws). 

+ In the clipping approximation, the diffraction loss is simply 2ir f^ drrp(r), where R is the mirror 
radius. In the ITM approximation, R is larger than all other length scales; however, the actual, finite 
value of R must be used in the clipping approximation for the diffraction loss to be nonvanishing. 
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Figure 4. The scaling of thermal noises with beam-width parameter D for mesa 
beams in the infinite-test-mass approximation. More specifically: a log-log plot of 
CiTu[r;p m esa.(D;r),p mesa ,(D ;r)] as a function of D/b. Here D /b = 1.76, which 
corresponds to a mirror radius of 12cm and 1 ppm diffraction losses. The curves are 
well approximated by power laws of the form C oc 1/ ' D 1 . 

4-1.3. Noise reduction by switching from a Gaussian beam to a mesa beam with the same 
diffraction loss and mirror radius Finally, the scaling law ( 1391) can be used to estimate 
the reduction in thermal noise by switching from a Gaussian beam to a mesa beam that 
has the same clipping-approximation diffraction loss on a mirror of the same radius. 

Two complications in the resized-beam scalings are not present when scaling from 
Gaussian to mesa beams. First, while the original and resized beams were associated 
with different-sized mirrors, now the Gaussian and mesa beams are associated with the 
same mirror. Second, when relating the Gaussian and mesa beams, there is no need to 
specify a fiducial beam size (i.e. there is no analogue of w Q and D Q ). Without these two 
complications, the Gaussian-to-mesa scaling is perhaps conceptually cleaner than the 
resized-beam scalings. 

Figure [5] shows Citm[t~', Pmesa(-D; r),p g!aiss (w, r )] for the sequence of beams shown in 
figure [2] (beams with 1 ppm diffraction loss in mirrors of the same radius R). The relative 
improvement in amplitude sensitivity increases monotonically with the mirror radius R, 
or equivalently, with D/w; however, when edge effects (i.e. finite-test-mass effects) are 
included, there is a limit to how much the sensitivity can be improved (section 14.2ft . 

4-2. Errors due to neglecting finite-test-mass (FTM) effects 

In the previous section, the ITM scaling laws predicted that, if the diffraction losses are 
held fixed, then the coating and substrate noises decrease monotonically with increasing 
beam width [figures [31 HI and [5]. In other words, for a given diffraction loss, the optimal 
beam width is simply "as large as possible." 

However, this conclusion is only as strong as the ITM approximation. Its validity 
can be checked by comparing the beam widths to the corresponding mirror dimensions. 
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Figure 5. The improvement in amplitude sensitivity when mesa beams 
are used instead of Gaussian beams. More specifically: a log-log plot of 
CnM[T;pmesa(D;r),p gaULSS (w]r)] as a function of D/w. For each mirror radius R, w 
and D are chosen so that the diffraction losses are 1 ppm in the clipping approximation. 
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Figure 6. How the Gaussian beam width parameter w compares to the mirror 
radius R and thickness H, when i) the radius R is fixed so the clipping-approximation 
diffraction loss is 1 ppm (unless a 10 ppm loss is indicated), and ii) the thickness H 
is then determined by holding the mass at 40 kg. the advanccd-LIGO baseline mirror 
mass. Each curve is proportional to w 1 . FS and Sap mean fused-silica and sapphire 
substrates. 



In our modelling, the mirror radii R are adjusted to maintain a constant clipping- 
approximation diffraction loss (CADL) [figure |2], while the thicknesses H is then 
determined by requiring the mirror mass be 40 kg — the design specification for Advanced 
LIGO. (Thus H will depend on whether the substrate is Fused Silica (FS) or sapphire 
(Sap), since the densities of these materials differ by a factor of about 2.) 

As shown in figures [6] and for the sequences of beam widths considered in 
section 14.11 w and D can approach or even exceed H while simultaneously being 
significant fractions of the R. Consequently, edge effects (finite test-mass effects) may 
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Figure 7. How the mesa beam width parameter D compares to the mirror radius R 
and thickness H, when i) the radius R is fixed so the clipping-approximation diffraction 
loss is 1 ppm (unless a 10 ppm loss is indicated), and ii) the thickness H is then 
determined by holding the mass at 40 kg. The mirror radius R for 1 ppm losses is 
shown on the top axis; the 10 ppm mirror radii are (from left to right) i? 10 PP m =13.94 
cm, 15.7 cm, 16.37 cm, 18.85 cm, and 21.36 cm. FS and Sap mean fused-silica and 
sapphire substrates. 



significantly change the sensitivity scalings depicted in figures El HI and [5j 

To estimate the importance of these edge effects, I compare the results in sections 
I4XT1 - 13X31 to the finite-test-mass (FTM) resultsH of Agresti and DeSalvo p}] (all 
types of thermal noise, 1 ppm CADL) and O'Shaughnessy, Strigin, and Vyatchanin [8] 
(substrate thermoelastic noise only, 10 ppm CADL). Specifically, from these data I read 
off the ratio 



/gFTM 

C F TM[r;Pu(r),Pk(r)} = J -^j. (46) 

This change in sensitivity can be compared to CiTM[r;pk(r),p u (r)}, the change in 
sensitivity obtained by the ITM approximation. Specifically, if 

^(rMlE ^f''"}! (47) 
CiTM[T;Pu{r),Pk{r)\ 

then |1 — A | is the fractional error made by using the ITM approximation to compute 
C[T;p u (r),p k (r)}. 

In the following subsections, I consider the errors |1 — A| made [sections 14.1.11 - 
14X3] by neglecting FTM effects. 



* The FTM data used here assume that the coating extends all the way to the edge of the substrate 
face. In advanced-LIGO, the coating radius will actually be several mm smaller than the substrate 
radius (the baseline substrate radius for advanced LIGO is 170 mm). 
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Figure 8. A log-log plot of A[T;p ga , uss {w; r),p gauss (w ; r)]. The fractional error of the 
sensitivity change made by neglecting edge effects is |1 — A|. Here w /b = 1.24, which 
corresponds to R = 12cm and 1 ppm diffraction losses. The FTM values are obtained 
by taking ratios of the noises calculated by Agresti and DeSalvo jTQ] ■ FS and Sap 
mean fused-silica and sapphire substrates. 



4-2.1. Resized Gaussian beam Figure M plots A[r; p gauss (u>; r),p gauss (u> ; r)] for mirror 
substrates made of fused silica, the baseline material for advanced LIGO mirrors [1] . For 
comparison, the figure also shows the corresponding values of A for sapphire substrates. 

When the substrate is fused silica, the ITM and FTM scaling laws agree to better 
than about 10% so long as R < 17cm, the advanced-LIGO baseline mirror radius [lj. 
As R increases beyond about 17 cm, |1 — A| increases dramatically (to about 50% when 
R = 21cm), because for such large radii the noise increases (e.g. [10l [8]) with R, while 
the ITM scaling laws predict [figure [3] that the noise always decreases with increasing 



When the substrate is sapphire, the FTM effects for the thermoelastic noises lead 
to errors that are comparable to the fused-silica FTM errors. For a mirror radius o^ 
R = 16cm, the fractional error 1 1 — A| for sapphire substrates is about 15% for substrate 
thermoelastic noise and about 20% for coating thermoelastic noise. 

4-2.2. Resized mesa beam The FTM effects in the resized-mesa-beam case are similar 
to the resized-Gaussian-beam FTM effects. Figure M plots A[r; p mGsa (D; r),p mesa (Z) ; r)]. 
When the substrate is fused silica and R < 17cm, the ITM scaling law errs by less than 
about 10% for the coating noises and by less than about 25% for the substrate Brownian 
noise. (The substrate thermoelastic noise is negligible when the substrate is fused silica 
[TO].) Again, the ITM scaling law disagrees more and more strongly as R is increased 
beyond 17cm. In this regime, the noise increases with R, but the ITM scaling law [figure 
[4] predicts that the noise always decreases with increasing R. 

jj When sapphire was the baseline test-mass material for advanced LIGO (it has since been abandoned 
in favour of fused silica), the baseline mirror radius was R = 15.7cm |21j . 



R. 
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Figure 9. A log-log plot of A[r;p mesa (D; r),p mesa (D ; r)]. The fractional error of the 
sensitivity change made by neglecting edge effects is 1 — A|. Here the diffraction losses 
are 1 ppm (unless 10 ppm is indicated), and D Q /b = 1.76 (D^pp 111 = 3.00), which 
corresponds to a mirror radius R — 12cm (i? 10 PP m = 13.94cm). The corresponding 
mirror radii are given on the top axis (1 ppm losses) and in figure [7] (10 ppm losses). 
The FTM values are obtained by taking ratios of the noises calculated by Agrcsti 
and DeSalvo [10], except for the 10 ppm values due to O'Shaughnessy, Strigin, and 
Vyatchanin [§]. FS and Sap mean fuscd-silica and sapphire substrates. (The fused- 
silica substrate thcrmoclastic noise is negligible; this case is omitted from the figure.) 
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When the substrate is sapphire, the FTM effects for the thermoelastic noises are 
comparable to the Brownian-substrate errors for fused silica. When R = 16cm, the FTM 
effects on the sapphire thermoelastic noises correspond to a fractional error |1 — A| of 
20% - 30%. 

4-2.3. Switching from a Gaussian beam to a mesa beam with the same diffraction loss 
and mirror radius The errors due to neglecting FTM effects in the Gaussian-to-mesa 
case behave qualitatively differently from (and are generally smaller than) the resized- 
beam errors. Figure [T0l plots A[r; p mcsa (-D; r),p gauss (w; r)] for fused silica and sapphire 
substrates. For both fused-silica and sapphire substrates, the coating sensitivity changes 
are not strongly sensitive to edge effects; in these cases, Cftm and Citm differ by less than 
about 10% even when the beam widths exceed 17 cm (and thus are significant fractions 
of R and H [c.f. figures |6] and [7j). The substrate sensitivity changes are more sensitive 
to edge effects, but even then the edge effects remain below about 15%, provided that 
R < 17cm for fused-silica substrates and R < 16cm for sapphire substrates. 

5. Conclusion 

Changing the shape of the laser beam in advanced LIGO can reduce the thermal noise, 
which is the limiting noise source at frequencies from 40 Hz to 200 Hz. In the Fourier 
domain, the relations between the thermal noise and the beam shape for semi-infinite 
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Figure 10. A log-log plot of A[r;p mosa (L>; r),p gauss (w : r)]. The beam width 
parameters w and D are chosen so that the diffraction loss is 1 ppm (unless 10 ppm 
is indicated). The corresponding mirror width for 1 ppm diffraction losses is shown 
on the top axis; the 10 ppm point corresponds to a mirror radius of 15.7 cm. The 
fractional error of the sensitivity change made by neglecting edge effects is |1 — A|. 
The FTM values are obtained by taking ratios of the noises calculated by Agresti and 
DeSalvo [10] . except for the 10 ppm value, which is due to O'Shaughnessy, Strigin, 
and Vyatchanin [8]. FS and Sap mean fused-silica and sapphire substrates. (The 
fused-silica substrate thermoelastic noise is negligible; this case is omitted from the 
figure.) 



mirrors take the form of simple scaling laws. Moreover, the coating thermal noises obey 
the same local scaling law. These results enable a straightforward comparison of the 
thermal noises for two different beam shapes when edge effects are neglected. The scaling 
laws predict the improvement of mesa-beam sensitivities vs. Gaussian-beam sensitivities 
quite well. For 40 kg, fused-silica mirrors, the substrate-noise scaling laws agree with 
the finite-mirror results within approximately 15% for mirror sizes not larger than the 
advanced-LIGO baseline size of about 17 cm; the coating-noise scaling laws agree with 
the finite-mirror predictions to better than about 10%. Therefore, the infinite-test-mass 
scaling laws may be a very useful tool for estimating optimal beam shapes for advanced 
LIGO and other future gravitational-wave interferometers. 
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In this appendix, I derive equation (I3T1) . which I use in the derivation of the scaling law 
(136]) for Brownian substrate noise. Then, I deduce equation (fT8j) . which I use in the 
derivation of the scaling law ( |23l) for Brownian coating noise. 
First, consider the integral 



drr [S 2 W - S w (6 - S z 



(A.1) 



Combining equations (113aj) and (11 3 el) gives 

e-S zz = ^- ^ 'dke- kz p(k) 
2/x Jo 



A + ft 



+ kz 



kJ (kr). 



(A.2) 



Inserting equations flA.2j) and (11 3d) into the left hand side of equation (1A. 1 1) yields 

drr [S 2 W - S w (0 - S, 

1 



V Jo 
where 



dk / dA;'e- (fc+fc ' )2 



-/'■ 



dr 



A + fi 
J^k^J^k'r) 



+ kz 



-/'■ 



A + /i 



P(k)P(k')I, 



k' 



dr Ji(kr) J (k V) . 



(A.3) 



(A.4) 



Since /c and are variables of integration, and since aside from I itself, ( 1A.3|) is 
unchanged by letting k *-> k', I can be rewritten as 



dr 



J 1 (kr)J 1 {k'r) 1„ 



o 



-fc' / dr J^Ax) Jo(ifc'r) - -k / dr Ji(fcV) J (fcr). (A.5) 



o 



o 



The integrals in (|A.5[) are special cases of equations (11.4.33), (11.4.34), and (11.4.42) 
of [22]: 

dr = — 77(A; - fc ) + ~^jV(k - k), 



drJi(kr)J (k'r) 



2k 
V {k-k') 



Here 77 is the unit step function. Inserting fll.6oj) and (11.6 £>[) into (IA.5I) shows that 







drr [5^ - S w (9 - S z 



0. 



which is equation (}3"Tj) . 

Next, combining equations (I16al) - fll6dj) shows that 



ocoat Q 

D W ~~ D W 1 2=0 
ncoat ocoat /, 



QCO&l ^OC 

Thus, setting z = in (1 1.7ft gives 



drr 

which is equation (1T81) . 



/ ocoat \ ^ ocoat Z/icoat ocoat \ 



;i.6a) 
(1.65) 

(1.7) 

;i.8a) 
(1.86) 

(1.9) 
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Appendix B. Junction conditions for the stress and strain of a statically 
deformed, semi-infinite mirror with thin coating 

The junction conditions (\16a\) - (I16dj) are listed in equation (A4) of [I] along with 
the statement that for these conditions to hold, the Poisson ratios of the coating 
and substrate should not be "too different." This restriction is actually unnecessary, 
provided that the coating is sufficiently thin. One can see this as follows: 

Because the coating adheres to the substrate surface, the substrate surface and 
coating have the same tangential displacement. Continuity of u r and immediately 
implies continuity of S rr and S w . A straightforward pillbox integration of the 
equilibrium condition VjTij = then shows that T zz and T rz are also continuous across 
the junction. 

All of the other junction conditions given in equation (A. 4) of [I] then follow, with 
one exception: the junction condition on S( rz ^ should read /i coat S'^ t = [J> snh S^, not 
S(°z) = S(rz)- But since T rz = on the coating surface (and thus also to high accuracy 
throughout the thin coating), this error is irrelevant; the correct junction condition is 
simply = = 0. 
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